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基于 DWT 最 优 多 子 图 和 SIFT 几何 校正 的 鲁 棒 水 印 算法 
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摘 要 : 为 提高 水 印 系 统 鲁 棒 性 ， 提 出 一 种 基于 PSO 算法 优化 的 DWT 多 子 图 水 印 算法 。 算 法 选择 DWT 变换 后 的 多 
个 子 图 作 DCT 变换 ， 然 后 利用 PSO 算法 寻找 最 优 的 多 子 图 组 合 权 重 ， 利 用 扩展 变换 QIM 原理 诺 入 和 提取 水 印 。 为 提 
高 算法 抵抗 几何 攻击 的 性 能 ， 提 出 基于 改进 SIFT 的 水 印 图 像 抗 几何 攻击 方法 ， 采 用 一 种 基于 圆 的 特征 点 描述 方法 改 
进 SIFT， 然 后 在 水 印 提 取 前 通过 SIFT 特征 点 匹配 进行 几何 校正 预 处 理 。 实 验 仿真 结果 表明 ， 在 保证 水 印 不 可 察觉 的 
基础 上 , 基于 PSO 算法 优化 的 DWT 多 子 图 水 印 算法 的 便 棒 性 优 于 基于 DWT 单一 子 图 的 水 印 算 法 ; 加 入 改进 SIFT # 
法 的 水 印 系 统 抵 抗 几 何 攻击 的 能 力 明显 提高 。 
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and SIFT geometric correction 
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Abstract: In order to improve the robustness of watermarking system, this paper proposed a DWT multi-subgraph based 
watermarking algorithm with PSO algorithm. The algorithm selected multiple subgraphs after DWT transform for DCT 
transform, then used PSO algorithm to find the optimal combination weights of multiple subgraphs. The algorithm embedded 
and extracted the watermark by using the spread transform QIM principle. To improve the performance of watermarking 
algorithm in resisting the geometric attacks, this paper proposed an anti-geometric attack method based on improved SIFT. It 
adopted a round based feature point description method to improve the SIFT algorithm. Before watermark extraction, it 
performed geometric correction processing by SIFT feature point matching. Experimental results show, based on the 
imperceptibility of watermark, the robustness of DWT multi-subgraph based watermarking algorithm with PSO algorithm is 
better than watermarking algorithm based on DWT single-subgraph. Watermarking system with improved SIFT algorithm can 
significantly improve the resistance to geometric attacks. 
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transform (SIFT) 
0 az 速 发 展 ， 其 为 数字 产权 保护 提供 了 一 种 安全 有 效 的 技术 手段 ， 
并 且 被 广泛 应 用 于 数字 媒体 领域 。 
21 世纪 以 来 ， 图 像 、 文 本 等 作品 的 数字 化 使 用 与 传播 越 来 根据 目前 国内 外 相关 研究 结果 ， 数 字 水 印 主要 分 为 变换 域 
越 普 及 ， 这 在 给 入 们 生活 带 来 便利 的 同时 也 带 来 了 隐患 一 电子 。 数字 水 印 和 空 (时 ) 域 数字 水 印 两 大 类 。 基 于 空间 域 的 水 印 算法 
盗版 越 来 越 猩 狐 。 数 字 水 印 技术 在 此 背景 下 应 运 而 生 并 得 以 迅 。” 直接 将 水 印信 息肉 入 到 载体 图 像 的 像素 信息 中 ， 变 换 域 算 法 则 
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录用 稿 李 莹 莹 ， 等 : 基于 DWT 最 优 多 子 图 和 SIFT 几何 校正 的 鲁 棒 水 印 算 法 


是 先 将 图 像 变 换 到 相应 的 变换 域 ， 再 在 变换 域 中 进行 水 印信 息 。 本 文中 提出 基于 改进 SIFT 的 水 印 图 像 抗 几何 攻击 算法 ,将 SIFT 
的 嵌入 。 基 于 变换 域 的 数字 水 印 技术 又 分 为 离散 傅 里 叶 变换 。 算法 中 基于 和 矩形 域 或 圆 环 的 特征 描述 方法 .3 改进 为 基于 圆 的 
(DFT)、 离 散 余 弦 变 换 (DCT) 和 离散 小 波 变换 (DWT)。 小 波 变换 。 ”特征 描述 方法 ， 在 保证 算法 几何 校正 准确 性 的 基础 上 大 大 节省 
是 一 种 时 频 的 局 域 变换 , 可 以 使 得 信号 分 析 在 不 同 尺度 上 进行 。 ”运行 时 间 ， 然 后 将 其 用 于 水 印 算法 中 来 提高 算法 的 抗 几 何 攻 击 
图 像 在 经 过 一 级 小 波 变换 后 , 会 得 到 四 幅 子 图 像 (图 1), 能力 。 


别 为 低频 子 图 LL 和 高 频 子 图 HH、LH、HL。 低频 子 图 中 保存 、 me 

了 原始 图 像 的 大 部 分 基本 特征 ， 多 级 小 波 分 解 后 得 到 的 最 高 层 | BWPDCTIRO SS 

的 低频 子 图 会 汇集 图 像 的 绝 大 部 分 重要 信息 ， 它 是 对 图 像 的 主 本 文 算法 在 仿真 过 程 中 使 用 的 载体 图 像 为 M*N 的 灰 度 图 
体 特 征 的 描绘 ， 因 此 这 部 分 也 被 称 为 逼近 子 图 ， 人 眼 视 觉 对 它 B KEJ m*n 的 二 值 图 像 (M28m,N>8n)。 基 于 量化 索引 调制 
相对 的 会 更 加 敏感 ， 其 他 三 个 高 频 子 图 保存 了 图 像 的 一 些 细节 ”(QIM) 的 水 印 算法 外 具有 较 好 的 鲁 棒 性 , 其 嵌入 信息 量 大 而 且 可 
zi A, , 是 对 图 像 边缘 信息 的 描绘 ,人 眼 对 这 部 分 的 敏感 度 略 低 ， 实现 盲 检测 ， 基 于 扩展 变换 的 量化 索引 调制 (ST-QIM) 算 法 将 扩 


但 是 因为 其 包含 图 像 的 边缘 信息 ， 高 频 部 分 发 生 的 细小 改动 又 思想 用 于 QIM 算法 ， 不 仅 具 有 QIM 算法 特点 ， 还 可 将 原本 
容易 影响 到 图 像 的 视觉 效果 中 3。 离散 余弦 变换 (DCT) 是 一 种 实 导 中 于 单一 系数 上 的 失真 进行 扩展 ， 具 有 更 好 的 安全 性 和 和 鲁 棒 
数 正 交 变换 ， 能 用 来 进行 去 相关 和 能 量 压缩 操作 ， 它 是 对 语音 。 TE. Li 和 Jiag 等 人 090 对 算法 中 的 量化 步 长 、 投 影 向 量 和 括 
和 图 形 信号 的 一 种 准 最 佳 变 换 。 基 于 DCT 的 数字 水 印 算法 计 ”入 点 选择 等 问题 进行 了 详细 的 研究 。 本 文 算法 将 利用 基于 扩展 
算 简单 且 易 于 实现 ,其 易于 将 水 印信 号 分 布 至 时 域 的 所 有 像素 ， ”变换 的 量化 索引 调制 算法 (ST-QIM) 实 现 水 印 杠 入 。 为 便于 下 文 
并 且 易 与 人 眼 视觉 模型 相 结 合 ， 兼 容 于 国际 数据 压缩 标准 ， 被 DR, HÆF DWT 单一 子 图 的 DWT-DCT 联合 变换 水 印 算法 


Eh 


iun 


研究 界 普遍 关注 山 , 结合 DCT 变换 的 解 相关 与 能 量 集中 性 能 和 II STOIM-DWCT-single. 
DWT 变换 的 时 频 与 多 分 辨 率 分 析 特 性 ，Cheng 等 人 B] 提 出 在 STOIM-DWCT.single 算法 结合 DCT 变换 与 DWT 变换 来 
DWTDCT 的 混合 域 中 嵌入 水 印 ， 并 用 logistics 序列 进行 水 印 ”嵌入 水 印 。 首 先 将 载体 图 像 进行 二 级 小 波 分 解 ( 图 1)， 取 二 层 
加 密 提高 系统 的 安全 性 。Saini 等 人 中 对 DWTDCT， 包 括 SVD 

ee da 中 单一 子 图 用 于 水 印 嵌入 ; 然后 将 该 区 域 进行 如 x 之 的 分 所 
分 解 的 各 种 组 合 域 水 印 柑 入 算法 进行 攻击 检测 仿真 分 析 。Singh 4m 4n 


等 人 提出 男 一 种 DWT-SVD 5 DCT 相 结合 的 水 印 算法 , 对 载 作 DCT 变换 。 
体 图 像 DWT 变换 后 的 水 平和 垂直 分 量 作 SVD 分 解 ， 水 印 图 
像 作 DCT 变换 后 利用 猫 映 射 置 乱 加 密 ， 然 后 将 加 密 水 印 骨 入 
载体 的 DWT-SVD 变换 系数 中 ， 提 高 水 印 系统 的 鲁 棒 性 与 安全 
性 。 


本 文 综合 DWT 和 DCT 变换 的 特点 ， 以 扩展 变换 QIM 为 
基础 算法 , 研究 基于 DWT 和 DCT 联合 变换 的 鲁 棒 性 数字 水 印 
技术 。 与 通常 的 DWT 或 DWTDCT 混合 域 水 印 算法 利用 某 个 


或 某 两 个 DWT 子 图 嵌入 水 印 的 方法 不 同 ， 本 文 在 研究 比较 基 图 1 小 波 分 解 图 

T DWT 不 同 子 图 的 水 印 算法 的 性 能 后 ， 还 提出 基于 DWT 多 水 印 嵌 入 流程 见 图 2。 有 具体 过 程 如 下 : 

子 图 的 水 印 算法 。 选 择 DWT 变换 后 的 多 个 子 图 作 DCT 变换 ， a) 分 别 读 取 载体 图 像 了 和 水 印 图 像 W. 

通过 一 定 的 权重 系数 将 其 组 合 起 来 用 于 水 印 柑 入 。 为 寻找 较 优 b) 采 用 Haar 小 波 将 载体 图 像 I 进 行 二 级 小 波 分 解 ， 得 到 多 
的 多 子 图 组 合 方法 ， 引 入 了 粒子 群 优 化 算法 (PSO) 四 ,设置 合理 ”个 分 辨 率 下 的 逼近 子 图 和 细节 子 图 ， 取 第 二 层 分 解 的 单一 子 区 


的 优化 函数 ， 通 过 粒子 群 优化 算法 的 迭代 运算 来 寻找 多 子 图 的  (LL2,HL2,LH2,HH2). 
最 优 组 合 。 中 将 对 应 子 图 进行 分 块 作 DCT 变换 ， 读 取 DCT 块 中 除了 

针对 水 印 的 攻击 , 常见 的 有 两 种 ;一 是 常规 信号 处 理 攻 击 ， (1,1) 的 其 他 系数 作为 向 量 x, (1,1) 系 数 用 来 决定 量化 步 长 , 然后 
即 噪声 攻击 、 滤波 攻击 、 压 缩 攻击 等 ; 二 是 几何 攻击 , 如 旋转 、 ”使 用 logistics 混沌 系统 来 产生 投影 向 量 v， 由 式 (1) ERR 
剪 切 、 尺 度 缩放 等 处 理 。 实 验 仿真 表明 提出 的 基于 粒子 群 优化 入 系数 y: 
的 DWT 多 子 图 水 印 算法 在 抵抗 常规 信号 攻击 时 表现 出 了 较 好 


c 


y-xv a) 
的 鲁 棒 性 ， 但 在 抵抗 几何 攻击 上 存在 不 足 ， 因 此 本 文 又 提出 加 
入 几何 校正 算法 对 受 攻击 载体 进行 预 处 理 .Lower7 提 出 的 SIFT qd) 利用 以 下 公式 将 水 印信 息 贬 入 子 图 的 系数 中 : 
算法 是 计算 机 视觉 中 一 种 经 典 的 匹配 算法 ， 具 有 比例 尺度 和 施 
转 不 变性 等 特点 ， 进 行 图 像 匹配 的 可 靠 性 非常 高 ， 但 该 算法 在 TUM €: 
性 能 方面 ， 如 计算 时 间 复 杂 度 、 匹 配 速 度 等 还 需 进一步 提升 。 me = round( qe) 68) 


:201804.02149v1 


chinaXiv 


录用 稿 
Oe = me — qe (4) 
if mod( me ,2) = w; 
E (5) 
Ju = me “A 


elseif mod(me,2) + w, 
f 6e>0 y,-(me-1).A (6) 
f oe<0 y,=(me+1):-A 
其 中 : ?为 待 嵌 入 系数 ;Jrw 为 帜 入 水 印 后 的 系数 ，A 为 步 长 ;wi 
为 水 印信 息 。 
e) 然 后 由 式 (7) 得 到 嵌入 水 印 后 的 DCT 系数 值 xw， 对 其 进 
fT DCT 反 变 换 得 到 小 波 系数 值 ， 再 对 得 到 的 系数 矩阵 进行 二 
级 小 波 重 构 ， 最 后 得 到 了 含水 印 的 载体 图 像 Iw 
x,-(y-y))vtx (7) 
根据 下 式 计算 峰值 信 品 比 : 


5 
PSNR-10 TES (8) 
MSE- v ON 1,6, DT (9) 


其 中 : ADI dj) 2) 380 ERAK E BU Jes RAE PUER BU RR EL s 


4 
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ChinaXi [A HN 


: 基于 DWT 最 优 多 子 图 和 SIFT 几何 校正 棒 水 


"yn 
Tee, 


将 每 个 二 层 子 图 均 做 分 块 DCT 变换 ， 就 形成 了 四 个 DCT 系数 
矩阵， 从 四 个 和 矩阵 中 选取 相应 的 DCT 系数 , 然后 对 每 个 系数 给 
予 一 个 权重 a1，a2，03 和 oa4， 将 其 组 合 起 来 ， 即 
y-a,DCT *a,: DCT,*a,- DCT, Ta DCT, (14) 

得 到 的 系数 y 即 为 待 量化 的 值 用 来 嵌入 水 印 。 用 扩展 变换 

的 思想 来 说 ， 基 于 DWT 单子 图 的 算法 中 的 扩展 变换 系数 来 
同一 个 DWT 子 图 的 DCT 变换 ， 而 基于 DWT 多 子 图 的 算法 中 
的 扩展 变换 系数 来 自 于 DWT 的 多 个 不 同 子 图 的 DCT 变换 ， 即 


y=x"v (15) 
其 中 : x S[DCT,, DCT ,, DCT}, DCT ,], DCT RBS i ^NT 
图 的 DCT 系数 值 ; y=[Q,0,,03,04]。 权重 ai, C2， aa 和 as 的 


选择 对 算法 的 性 能 至 关 重 要 ， 本 文 将 通过 粒子 群 算法 来 搜索 其 
最 优 解 。 粒 子 群 优 化 算法 (particle swam optimization, PSO) 是 一 
种 全 局 优化 算法 ， 算 法 通过 模拟 鸟 类 飞行 砚 食 时 相互 之 间 的 群 
体 协作 行为 使 得 整体 最 优化 。 在 PSO 优化 算法 中 ， 待 优化 题目 
的 所 有 潜在 解 被 看 做 是 搜索 空间 内 的 一 群 鸟 ， 每 个 解 相 当 于 一 
只 鸟 ， 每 只 鸟 都 由 一 个 矢量 来 决定 它 飞 行 的 速度 与 方向 即 速 度 
矢量 ， 在 算法 中 解 也 就 是 粒子 。 根 据 实际 问题 确定 对 应 的 优化 


MSE 为 载体 图 像 典 入 水 印 前 后 的 均 方 误差 。 理 论 上 PSNR 值 在 
40 dB 以 上 ， 水 印 图 像 已 具有 较 好 的 不 可 见 性 。 


水印 的 提取 过 程 是 舱 入 的 道 过 程 。 
PA 
dp = " (10) 
m, — round(q, ) (11) 
w,'= mod(m, ,2) (12) 
比特 误 码 率 (bit errorratio,BER ) 用 来 衡量 提取 水 印 的 恢复 
情况 ， 其 值 越 小 说 明 水 印 恢复 的 越 好 ， 计 算 方 法 为 
BER= 传 输 中 的 误 码 /传输 的 总 码 数 (13) 


二 级 离散 
体 图 一 | DCT 变换 
小 波 变换 e| 


嵌入 水 印 后 le 二 级 小 a DCT 反 变换 
的 载体 图 像 波 重 构 


图 2 基于 DWT 单一 子 图 的 DWT-DCT 联合 变换 水 印 算法 

为 便于 区 分 ,将 基于 二 次 DWT 分 解 后 LL2 子 图 的 STQIM- 
DWCT-single 算法 ， 也 是 文献 [3] 中 使 用 的 联合 变换 方法 ， 记 做 
DWCT(LL2)， 采 用 其 他 三 个 高 频 子 图 的 STQIM- DWCT-single 
45:4) 37g DWCT(HL2), DWCT(LH2), DWCT(HH2). 


ra 


2 “基于 粒子 群 算法 优化 的 DWT 多 子 图 水 印 算法 
与 之 前 的 基于 DWT 单一 子 图 的 联合 变换 水 印 算法 不 同 ， 

本 文 提出 选择 DWT 的 多 个 子 图 用 来 水 印 嵌 入 。 本 文 提 出 选择 

慨 逼 近 子 图 进行 再 分 解 后 得 到 的 二 层 四 个 子 图 用 来 水 印 嵌 入 ， 


函数 来 决定 适应 值 ， 根 据 适 应 度 值 寻找 最 优 粒子 。 粒 子 通 过 跟 
踪 局 部 和 全 局 的 最 优 粒子 更 新 自己 ， 通 过 迭代 的 方法 找到 最 终 
的 最 优 解 外 。 粒 子 群 算法 的 运算 规则 比 遗 传 算法 简单 ， 它 不 同 
于 遗传 算法 中 的 “交叉 ”和 “变异 ”， 而 是 通过 已 搜索 到 的 最 
优 值 来 寻求 全 局 最 优 解 。PSO 算法 因 其 算法 简单 易于 实现 、 收 
敛 精度 高 等 得 到 了 研究 界 的 广泛 关注 。 

PSO 算法 中 几 个 重要 定义 : 粒子 I ENA 
X pi—(prs pis v py) 飞行 速度 矢量 为 si— Ut SN)o 
每 个 粒子 根据 目标 函数 计算 相应 的 适应 度 值 Gimess), 根据 适 
度 值 记录 单个 粒子 的 当前 位 置 p; 和 到 目前 为 止 所 发 现 的 最 好 位 
置 (pbesti))。 此 外 ， 还 要 记录 群体 粒子 中 截止 目前 所 找到 的 全 局 
最 优 位 置 (gbest)。 然 后 粒子 通过 自己 单个 的 最 好 记录 和 群体 的 
最 好 记录 来 决定 下 一 步 的 进化 。 

PSO 算法 是 一 种 迭代 算法 ， 它 需要 先 初始 化 一 群 随机 解 作 
为 粒子 群 ， 然 后 进行 迭代 运算 来 寻找 到 最 优 解 。 在 每 一 次 迭代 
中 ， 粒 子 通 过 与 两 个 最 好 位 置 (pbest;，gbest) 比 较 来 更 新 自己 。 


空间 的 位 置 矢量 


(si, 52; 


k 


粒子 通过 下 面 的 两 个 公式 来 进行 速度 和 位 置 的 更 新 : 
S; 2 5; +c: rand() : (pbest; — p;) +c, :rand(). (gbest— p;) (16) 
pi=Pits, i-lA..N (7) 


其 中 : % 是 粒子 的 速度 矢量 ; 产 代 表 粒 子 当 前 的 位 置 ; pbest; 和 
gbest 如 上 述 定义 ; cl 和 co 是 学 习 因 子 。 对 于 每 个 粒子 ,在 每 一 
维 中 都 会 对 其 速度 设置 一 个 阔 值 V mas (Vina >0) 来 限制 粒子 的 速 
度 ， 使 得 粒子 的 速度 不 超过 最 大 速度 〈 或 最 小 速度 ) 。 

为 便于 计算 ,将 基于 PSO 优化 的 DWT 多 子 图 的 DWT-DCT 
联合 变换 水 印 算 法 记 为 STQIM-DWCT-PSO multi。STQIM- 
DWCT-PSO multi 算法 流程 框图 见 图 3。 
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过 滤 掉 不 稳定 的 极 值 点 。 描 述 特征 点 时 ， 首 先 计算 特征 点 的 像 
RABB | 2 pwr ER 素 精度 位 置 ， 再 在 高 斯 图 像 上 统计 特征 点 特定 邻 域内 像素 的 梯 


"rum 度 方向 与 梯度 和 ， 确 定 特征 点 的 主 方向 。 描 述 特征 点 时 ， 利 用 

关键 点 所 在 的 尺度 空间 (高 斯 金字 塔 的 某 一 层 ) 确定 特征 点 周 

[— > kmga T ] 围 的 16 个 4*4 UD AUR, 计算 每 个 子 域 的 梯度 方向 直方 图 ( 梯 

EPRE Y ds =e 度 方向 均匀 分 为 m 个 方向 ,通常 m 29 80 ; 然后 将 计算 出 的 方 

i mn Y Y 向 梯度 直方 图 按 位 置 进行 排序 ， 进 而 构成 了 16*8 即 128 维 的 

^ 水 印 提取 | | 水 印 提取 || Amm 特征 描述 子 向 量 ， 最 后 对 特征 描述 子 向 量 进行 归 一 化 处 理 ， 以 

a "S me E 减少 光照 亮度 变化 的 影响 。 图 像 匹 配 时 通常 通过 计算 特征 点 之 
闻 的 欧 氏 距离 来 进行 图 像 特 征 点 的 匹配 判定 。 


V 
sss |E 经 典 SIFT 算法 保持 有 较 高 的 匹配 准确 率 ， 但 因为 要 生成 
128 维 的 描述 子 ， 不 但 要 耗费 较 多 的 时 间 ， 还 占用 较 大 的 内 在 


空间 ， 当 特征 点 较 多 时 严重 影响 特征 匹配 的 速度 。 之 后 研究 者 
们 提出 了 多 种 SIFT 改进 算法 。 文 献 [8] 提 出 一 种 在 极 坐标 下 的 
图 3 基于 PSO 优化 的 水 印 算法 流程 特征 描述 方法 ， 以 特征 点 为 圆心 扩散 出 多 个 圆 ， 每 个 圆 拆 分 成 

PSO 算法 优化 水 印 系 统 的 过 程 如 下 : 四 个 扇形 ;然后 以 x 坐标 为 界 进行 垂直 镜面 反射 ， 得 到 新 的 特 

a) 初 始 化 m 个 系数 组 (al，a2，a3，o4) 作 为 粒子 群 初始 位 。” 征 描述 区 域 计算 特征 描述 子 。Amin 等 人 0 提出 一 种 自 适 应 分 
lL pi， 并 产生 一 组 随机 粒子 作为 变换 速度 sio 区 的 SIFT 改进 算法 ， 在 对 数 极 坐标 网 络 中 ， 划 分 三 个 圆 形 区 

b) 初 始 化 所 有 pbesti 和 gbest 的 适应 度 值 为 0。 或 ， 在 每 个 圆 形 区 域 中 划分 出 多 个 扇形 区 域 ， 在 每 个 扇形 区 域 

9 由 峰值 信 噪 比 和 误 码 率 设 计 适 应 度 函 数 ， 根 据 水 印 算法 ”内 统计 不 同 的 梯度 方向 直方 图 来 描述 特征 点 。 卢 官 明 等 人 DO] 提 
进行 水 印 嵌 入 与 提取 并 计算 每 个 粒子 的 适应 度 值 fitness。 出 基于 圆 环 域 的 特征 描述 改进 方法 ， 以 特征 点 为 圆 形 ， 向 外 扩 

dd) 每 个 粒子 都 将 其 当前 的 适应 度 值 与 本 粒子 对 应 的 最 好 位 。 散 出 多 个 圆 形 区 域 ， 计 算 每 个 圆 形 区 域内 的 梯度 直方 图 来 描述 
(pbesti) 的 适应 度 值 比较 ， 如 果 其 结果 更 好 ， 则 更 新 该 粒子 的 。” ”特征 点 。 但 它们 或 是 在 提高 描述 质量 时 未 考虑 算法 运行 时 间 ， 
最 好 位 置 pbesti 为 当前 位 置 。 或 是 在 提高 算法 速度 时 影响 了 匹配 的 准确 率 。 

6) 将 每 个 粒子 的 适应 度 值 与 整个 群体 粒子 经 过 的 最 好 位 置 本 文 对 SIFT 算法 的 特征 描述 子 提出 一 种 新 的 改进 ， 同 样 
(gbesb 的 适应 度 值 作 比较 ， 如 果 某 个 粒子 的 结果 较 好 ， 则 将 其 ”采用 同心 圆 的 形式 ,但 不 同 于 统计 圆 环 内 点 的 方向 梯度 直方 图 ， 
位 置 作 为 当前 整体 的 最 好 位 置 gbest。 而 是 选择 圆 上 的 像素 点 来 统计 其 方向 梯度 直方 图 。 改 进 算法 首 

根据 式 (16)(17) 调 整 粒 子 的 速度 和 位 置 。 先 以 特征 点 为 圆心 ， 选 取 4 个 不 同 半径 的 同心 贺 上 的 点 作为 采 

判断 是 否 达 到 结束 条 件 ， 没 有 达到 结束 条 件 则 跳 至 c); 样 点 ， 在 同一 个 贺 上 的 采样 点 之 间 上 共有 相同 的 间距 ， 每 个 贺 上 


若 满 足 条 件 ， 则 此 时 的 gbest 即 为 最 优 的 系数 组 。 的 点 在 特征 点 的 指定 方向 上 。 然 后 统计 相同 半径 圆 上 的 点 的 方 
常用 的 算法 迭代 终止 方法 有 两 种 ， 一 是 确定 最 大 友 代 次 数 向 梯度 直方 图 ， 根 据 圆 的 不 同 对 其 进行 排序 继而 构成 了 4*8 即 
G; 二 是 确定 最 小 (或 最 大 ) 适应 阔 值 ， 当 粒子 群 搜索 到 的 最 优 ， ”32 维 的 特征 描述 子 向 量 。 本 文 算法 选择 与 矩形 邻 域 一 样 个 数 的 
位 置 满足 预定 的 适应 阔 值 则 迭代 终止 。 本 文中 选择 最 大 友 代 数 ”64 个 采样 点 (图 4) ， 采 样 点 的 确定 方式 如 下 : 设 采 样 点 的 坐 
作为 终止 条 件 。 标 为 (XG),Y())771,2,..,64, dp; 


、 - r(i) - 4(14 ((i-1)mod 4)) (18) 
3 ”基于 改进 SIFT 的 水 印 图 像 抗 几何 攻击 方法 un 
w—— " H -- 9) = fix( ——) (19) 
为 了 提高 数字 水 印 算法 抗 几何 攻击 的 鲁 棒 性 ， 近 年 来 研究 4 
者 们 提出 了 许多 抵抗 几何 攻击 的 水 印 方案 ， 其 中 包括 几何 不 变 . f 27 (i) 
" PROP re 2 . i X (i) 2 r(i) cos( ———) (20) 
域 方 法 、 和 矩 估 计 方 法 和 基于 特征 点 的 方法 等 。 其 中 较为 典型 的 1 
是 基于 特征 点 的 SIFT 算法 。SIFT 算法 是 计算 机 视觉 中 一 种 有 , — . 2mp (i) 
: Y (i) = r(i) sin( ———) Q1) 
效 的 匹配 算法 ， 具 有 旋转 不 变性 和 比例 尺度 不 变性 等 特点 。 16 
SIFT 算法 是 局 部 特征 提取 算法 , 它 通过 在 不 同 尺 度 空间 内 计算 


极 值 点 ， 并 将 对 应 的 位 置 、 尺 度 和 主 方向 等 特征 量 提 取出 来 。 
SIFT 的 特征 点 是 通过 将 差分 高 斯 图 像 中 的 每 一 个 像素 点 与 当 
前 尺度 上 和 相 邻 的 上 、 下 层 尺 度 上 的 26 个 邻 域 点 相 比较 而 得 
到 的 极 大 值 或 极 小 值 ， 并 通过 泰勒 展开 式 和 Hessian 矩阵 计算 
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图 4 改进 的 特征 描述 子 采样 模式 


加 入 改进 SIFT 进行 几何 校 1 


见 图 5。 具 体 步骤 如 下 : 读 取 载 体 图 


E 预 处 理 的 水 印 算法 流程 框 


多 


合 变换 , 组合 多 个 DWT 子 图 的 DCT 系 


嵌入 ， 嵌 入 水 印 后 的 系数 先进 行 DCT 反 
含水 印 的 载体 图 像 。 对 藤 入 水 
与 描述 ， 保 存 这 些 特征 


IFT 特征 点 提取 


ST-QIM 算法 进行 水 印 
变换 再 进行 小 波 重 构 ， 最 后 得 到 
印 后 的 载体 图 像 进行 S 

点 作为 密 钥 用 于 后 期 校 

各 种 攻击 , 在 水 印 


形 了 对 其 进行 DWTDCT 


AE 


ZUR T RUNZKED, A 


dH 


正 。 经 过 传输 使 用 后 的 载体 图 像 会 受到 
提取 之 前 , 先 对 受 攻击 载体 进行 SIFT 特征 点 


提取 与 描述 ， 然 后 与 保存 的 特征 点 密 钥 进行 匹配 ， 判 断 图 像 是 


ChinaXiv 合 作 期 刊 
Tee, ar pwrR& S QE ELLE 法 


fitness = PSNR +u 23 (1- BER (i)) Q2) 


其 中 : / 为 参数 ; BER(i) 指 不 同 攻 击 后 提取 出 水 印 的 误 码 率 。 当 
PSNR 值 愈 大 ， 表 示 水 印 图 像 的 不 可 见 性 越 好 ; BER HAN, 
表示 恢复 出 的 水 印 图 像 质量 越 好 。 综 合 两 个 指标 ， 当 适应 度 值 
fitness 越 大 ， 表 示 算 法 综合 性 能 越 好 。 
为 验证 所 提 算 法 性 能 ， 本 文 将 基于 DWT 单一 子 图 的 四 个 
算法 , 即 DWCT(LL2) (文献 3]), DWCT(HL2), DWCT(LH2), 
DWCT(HH2)， 与 本 文 提出 的 基于 PSO 算法 优化 的 DWT €T 
图 水 印 算法 (DWCT-PSO) 进 行 比较 。 为 了 便于 比较 算法 鲁 棒 性 ， 
通过 调节 水 印 嵌入 强度 ,将 所 有 算法 的 峰值 信 品 比 调节 至 45 dB, 
以 保证 巾 入 水 印 图 像 具 有 相同 的 不 可 见 性 ， 然 后 对 其 鲁 棒 性 进 
行 比较 。 仿 真 结果 见 图 6。 
从 图 6 对 几 种 算法 进行 仿真 得 到 的 结果 对 比 可 以 看 出 ， 
T DWT 单个 子 图 的 四 种 算法 中 ， 选 择 逼 近 子 图 用 于 水 印 供 入 
的 算法 在 综合 鲁 棒 性 上 优 于 其 他 三 种 算法 。 而 本 文 提 出 的 基于 
DWT 多 个 子 图 , 并 经 过 PSO 算法 优化 的 DWCTPSO 算法 效果 


否 经 历 了 几何 攻击 ， 确 定 受到 几何 攻击 的 类 型 与 尺度 ; 然后 进 
行 相应 的 校正 处 理 ， 最 后 再 对 校正 后 的 载体 图 像 进 行 水 印 提 取 
操作 。 
载体 DWTDCT STQIM En DCT 反 变 换 
Er 联合 变换 水 印 嵌入 与 小 波 重 构 | 8KM 
载体 x 
水 名 图 人 SIFT 特征 点 
msi Xn 
Senis x" 
提取 出 STQIM DWTDCT SIFT 特征 点 
的 水 印 水 印 提取 联合 变换 提取 与 描述 
图 5 基于 改进 SIFT 的 抗 几何 攻击 水 印 算法 流程 


4 ”实验 仿真 及 分 析 


本 文 的 仿真 实验 中 选取 多 幅 分 辨 率 为 256*256 的 灰 度 图 像 
“ 东 学 ”字样 的 二 值 图 像 作为 水 印 图 像 。 


作为 载体 图 像 ， 
水 印 嵌 入 与 提取 算法 以 ST-QIM 算法 为 基础 。 


"UH 


41 基于 粒子 群 算法 优化 的 DWT 多 子 图 水 印 算法 


仿真 实验 时 ， 基 于 DWT 单子 图 的 DWT-DCT 联合 变换 水 


印 算法 中 ， 首 
二 层 子 图 中 的 


先 将 载体 


HT E 


FH 3 


对 其 仿真 结 
合 变换 水 印 算 
的 DCT 变换 ， 


果 进 行 比 
法 中 提出 
然后 结合 


的 子 图 组 合 系 


多 的 系数 因子 。 


图 像 进 行 二 级 DWT 变换 ， 然 后 只 选择 
进行 相应 的 DCT 变换 用 来 水 印 岩 入 ， 
较 。 基 于 DWT 多 子 图 的 DWT-DCT 联 


选择 二 层 子 图 中 的 多 幅 子 图 ， 均 作 相应 
多 个 子 图 用 来 水 印 杠 入 。 为 了 选择 合适 


数 ,首先 需要 由 粒子 群 优化 (PSO) 算 法 来 求解 出 最 


本 文中 利用 峰 


数 ， 即 


在 PSO 算法 中 ， 需 要 选择 合适 的 适应 度 函数 ， 


值 信 噪 比 (PSNR) 和 误 码 率 (BER) 来 构造 适应 度 函 


最 优 ， 其 不 仅 保证 水 印 的 不 可 见 性 ， 在 抵抗 各 种 信号 攻击 的 鲁 
棒 性 上 明显 优 于 基于 DWT 单个 子 图 的 水 印 算法 。 


Dworu 
+ WCG 
DWCT(LH2) 
o.s- DWCT(HH2) 
DWCT-PSO 


0.6|- 


ot = : : : : : : : 
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
variance of Gaussian noise x 10? 
r Herna = -+ 
(a) 高 斯 噪声 攻击 
0.8 E 
DWCT(LL2) 
07|- t— DWCT(HL2) 
DWCT(LH2) 
ER DWCT(HH2) 
E 一 二 一 DWCT-PSO 


Qi : T 5 eL 
20 30 40 50 60 70 80 90 100 
quality factor of JPEG 


(bJPEG 压缩 攻击 
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Be gp 表 2 特征 提取 时 间 比 较 


ab | DWCT(LL2) 
aal, i -— 特征 点 提取 时 间 /s 
oss. E 文献 [7] 文献 [13] 改进 SIFT 
i 尺度 缩放 0.5 1.9 1.3 0.5 
| 尺度 缩放 1.5 93 6.5 43 
尺度 缩放 2.0 21.2 17.5 5.4 
旋转 2? 7.0 4.6 1.2 
o MMC TN EE 旋转 15° 8.5 5.8 1.6 
apeppernoge o o CO 旋转 30。 Ms rs 30 
(c) 椒 盐 噪 声 攻击 旋转 60* 7.9 5.3 3.1 
图 6 不 同 算法 鲁 棒 性 比较 旋转 120" 11.1 7.6 3.1 
图 7 和 8 分 别 是 利用 DWCT(LL2) 算 法 即 文献 [3] 算 法 与 平移 (50,0) 8.9 5.9 1.4 
DWCT-PSO 算法 进行 水 印 葵 入 后 的 含水 印 载体 和 通过 提取 算 平移 (0,50) 5.8 3.7 2.4 
法 提取 出 的 水 印 图 像 。 平移 (50,50) 4.7 3.1 1.2 


Watermarked Image 


表 1 和 2 中 己 有 的 SIFT 算法 与 本 文 提出 的 改进 SIFT 算 
法 的 几何 攻击 校正 和 描述 子 生成 时 间 的 比较 可 以 看 出 ， 改 进 的 


Fr. SIFT 算法 仍 保持 较 高 的 校正 准确 性 , 但 描述 子 生 成 时 间 相 比 于 
e 文献 [7] 缩 得了 近 四 倍 。 这 是 因为 对 于 特征 点 的 描述 ， 越 靠近 特 


征 点 的 邻近 点 对 描述 的 贡献 越 大 ， 所 以 取 点 时 在 距离 特征 点 最 
近 的 圆 上 取 点 最 为 密集 ;为 了 增加 描述 准确 性 ， 将 原 算法 中 邻 


图 7 文献 [3] 算 法 (LL2 子 图 )， 仿 水印 载体 ( 左 )， 提 取水 印 ( 右 ) 近 域 由 16 扩大 为 32， 从 而 保证 了 较 高 的 校正 准确 性 。 因 为 只 
ee N 需 计 算 四 个 同心 圆 上 的 点 ， 且 每 个 圆 上 的 所 有 点 距离 特征 点 的 
方向 与 距离 都 是 一 定 的 ， 大 大 减少 了 运算 量 。 同 时 ， 每 个 特征 
Fr. -本 点 的 特征 描述 向 量 由 128 维 减 少 至 32 维 ， 也 大 大 节省 了 程序 

pe" 运行 所 需 的 内 存 空间 。 
] 将 提出 的 SIFT 改进 算法 与 水 印 算法 结合 使 用 ， 载 体 图 像 
"1E EP n 嵌入 水 印 后 ， 先 由 改进 SIFT 算法 提取 出 特征 点 并 记录 下 特征 
图 8 DWCTPSO 算法 ， 含 水 印 载体 £) ， 提 取水 印 〈 右 ) 描述 子 用 于 后 期 校正 。 在 水 印 提取 时 ， 算 法 首先 对 受 攻击 载体 
4.2 基于 改进 SIFT 的 水 印 图 像 抗 几何 攻击 方法 图 像 进行 特征 点 提取 ， 然 后 将 其 与 之 前 保存 的 特征 点 及 描述 子 
表 1 和 2 给 出 了 经 典 SIFT 算法 (文献 [7])、 文 献 [13] 中 的 基 ，” 作 匹配 ， 以 此 来 判断 含水 印 载体 图 像 是 否 受 到 几何 攻击 并 判断 


于 圆 环 域 的 改进 SIFT 算法 与 本 文 提出 的 改进 SIFT 算法 的 抗 几 ”出 几何 攻击 类 型 。 若 图 像 受到 了 几何 攻击 ， 则 先 对 含水 印 载体 


何 攻击 的 校正 检测 预 估 结果 与 特征 提取 时 间 的 比较 。 进行 几何 校正 ， 然 后 再 提取 水 印 。 图 9~11 为 未 引入 SIFT 校正 
表 1 几何 校正 预 估 结 果 比 较 的 水 印 算法 与 本 文 抗 几 何 攻击 水 印 算法 的 比较 。 左 一 为 受到 几 
几何 校正 预 估 结 果 何 攻 击 后 的 载体 图 像 左 二 为 未 进行 几何 校正 提取 出 的 水 印 图 像 。 
文献 [7] 文献 [13] 改进 SIFT 由 图 9~11 可 看 出 , 与 未 使 用 几何 校正 的 水 印 算法 相 比 , 基于 改 
尺度 缩放 0.5 0.4997 0.4995 0.4995 xt SIFT 的 抗 几何 攻击 水 印 算 法 能 较 准确 的 校正 图 像 并 提取 出 
尺度 缩放 1.5 1.4970 1.4981 1.4993 水 印 。 
尺度 缩放 2.0 2.0010 2.0014 2.0008 
旋转 2° 2.0710 2.0607 2.0462 
旋转 15° 15.014 14.995 14.992 
旋转 30* 30.033 30.052 30.013 
旋转 60* 60.017 60.063 60.022 
旋转 120。 119.998 120.016 120.008 图 9 图 像 放大 2 倍 ， 未 校正 ( 左 二 ), 图 像 较 正 及 提取 的 水 印 ( 右 一 二 ) 
平移 (50,0) (50,0.017) (49.9804,0) (50,0.018) 
平移 (0,50) (0,49.985) (0,49.985) (0.279,49.41) 


平移 (50,50) (50.01,49.99) (50.027,49.95) ^ (50.30,49.90) 
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图 10 图 像 平移 ,未 校正 ( 左 二 ), 图 像 较 正 及 提取 的 水 印 ( 右 一 二 ) 


50 — 100 150 200 250 


旋转 5 后 图像 


图 11 图 像 旋转 2”, 未 校正 ( 左 二 ) 图 像 较 正 及 提取 的 水 印 ( 右 一 二 ) 
同时 , 为 了 验证 本 文 提 出 的 抗 几何 攻击 水 印 算法 的 有 效 性 ， 
将 本 文 算 法 与 Fazli 等 人 0 提出 的 算法 进行 抗 几何 攻击 仿真 比 
较 。 为 进行 鲁 棒 性 比较 , 将 基础 水 印 算法 的 PSNR 值 设置 为 57 
dB. X 3 给 出 了 分 别 利 用 两 种 算法 进行 几何 校正 后 提取 出 的 水 


印 和 原水 印 的 NC 值 。 由 表 中 数据 可 看 出 ， 虽 然 本 文 算法 抗 几 
何 攻 击 的 鲁 棒 性 稍 逊 色 于 文献 [1 和 ,但 恢复 出 水 印 的 平均 NC 值 


在 0.95 以 上 。 通 常 利用 水 印 算法 恢复 出 的 水 印 NC 值 在 0.7 以 
上 ， 则 认为 该 水 印 算法 是 行 之 有 效 的 中 ， 因 此 本 文 提 出 的 基于 
改进 SIFT 的 抗 几 何 攻 击 水 印 算法 是 一 种 有 效 的 ,校正 精度 高 、 

性 能 稳定 的 抗 几何 攻击 水 印 算法 。 
表 3 本 文 算法 与 文献 [1 和 1 算法 的 抗 几何 攻击 性 能 比较 (NC) 


攻击 文献 [14] 本 文 算法 
尺度 缩放 0.5 0.980 0.965 
尺度 缩放 1.5 0.986 0.985 
旋转 307 0.975 0.936 
旋转 80° 0.977 0.973 
平移 (3525) 0.987 0.983 
中 心 剪 切 20% 1 0.98 
5 ”结束 语 


g! 


本 文 在 DWT 域 和 DCT 域 图 像 水 印 算法 基础 上 , 研究 基 了 
DWT 与 DCT 联合 变换 的 水 印 算法 。 实 验 仿真 结果 表明 ， 基 于 
DWT 单一 子 图 的 水 印 算法 中 ， 利 用 通 近 子 图 众 入 水 印 的 算法 
性 能 最 优 : 本 文 提 出 的 基于 粒子 群 算法 优化 的 DWT 多 子 图 水 
印 算法 的 鲁 棒 性 明显 优 于 基于 DWT 单一 子 图 的 水 印 算法 。 针 
对 所 提 算 法 在 抵抗 几何 失真 方面 的 缺陷 ,本 文 提出 利用 SIFT 算 
法 进行 几何 校正 预 处 理 , 并 针对 SIET 算法 的 不 足 进行 了 改进 。 
改进 后 的 SIFT 算法 用 于 接收 图 像 的 几何 校正 预 处 理 ， 不 仅 具 
有 较 高 的 校正 准确 率 ， 算 法 运行 时 间 也 大 大 减少 了 ， 还 降低 了 
程序 运行 所 需 的 内 存 空间 。 将 改进 的 SIFT 算法 应 用 于 水 印 算 
法 中 对 含水 印 载体 图 像 进行 几何 攻击 判断 并 校正 ， 明 显 提高 了 
水 印 提取 的 准确 率 ， 是 一 种 有 效 可 行 的 水 印 算法 。 综 上 所 述 ， 
本 文 所 提 的 水 印 算法 是 一 种 讶 检测 的 鲁 棒 性 数字 水 印 方法 ， 对 
水 印 有 较 好 的 不 可 见 性 ， 对 常见 的 攻击 具备 较 强 的 鲁 棒 性 ， 
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可 推广 至 彩色 图 像 水 印 算法 中 。 下 一 步 , 如 何 将 SIFT 算法 与 水 
印 算法 更 有 效 地 结合 为 一 体 ， 实 现 更 加 简洁 高 效 的 水 印 检测 系 
统 仍 需 要 进行 深入 的 研究 。 


参考 文献 : 


[1] E4, KR, EAA., 数字 水 印 原理 与 技术 [M]. 北京 : 科学 出 版 社 ， 
2007: 63-65. (Wang Ying, Xiao Jun, Wang Yunhong. Digital watermarking 
principles and techniques [M]. Beijing: Science Press, 2007: 63-65. ) 

[2] 4&4. 小 波 变换 与 图 像 处 理 [M] 北京 : 中 国 科学 技术 大 学 出 版 社 ， 
2010: 15-24. (Ni Lin. Wavelet Transform and Image Processing [M]. Beijing: 
University of Science and Technology of China Press, 2010: 15-24. ) 

[3] Cheng Wei, Li Zhaodan. Robust watermarking algorithm of color image 
based on DWT-DCT and chaotic system [C]// Proc of IEEE International 
Conference on Computer Communication and the Internet. 2016: 370-373. 

[4] Saini LK, Shrivastava V. Analysis ofattacks on hybrid DWT-DCT algorithm 
for digital image watermarking with MATLAB [J]. International Journal of 
Computer Science Trends and Technology, 2014, 2 (3): 123-126. 

[5] Singh D, Singh S K. DWT-SVD and DCT based robust and blind 
watermarking scheme for copyright protection [J]. Multimedia Tools and 
Applications, 2017, 76 (11): 13001-13024. 

[6] 思 格 尔 伯 里 特 . 计算 群体 智能 基础 [M]. 谭 营 , 译 . 北京 : 清华 大 学 出 
版 社 ，2009: 65-89. (Engelbrecht A P. Fundamentals of Computational 
Swarm Intelligence [M]. Tan Ying, Trans. Beijing: Tsinghua University 
Press, 2009: 65-89. ) 

[7] Lowe D G. Distinctive image features from scale-invariant keypoints [J]. 
International Journal of Computer Vision, 2004, 60 (2): 91-110. 

[8] Liao Kaiyang, Liu Guizhong, Hui Youshi. An improvement to the SIFT 
descriptor for image representation and matching [J]. Pattern Recognition 
Letters, 2013, 34 (11): 1211-1220. 

[9] Chen B, Wornell G W. Quantization index modulation: a class of provably 
good methods for digital watermarking and information embedding [J]. 
IEEE Trans on Information Theory, 2001, 47 (4): 1423-1443. 

[10] Li Xinchao, Liu Ju, Sun Jiande, et al. Step-projection-based spread 
transform dither modulation [J]. IET Journal on Information Security, 2011, 
5 (3): 170-180. 

[11] Jiang Yanling, Zhang Yifeng, Pei Wenjiang, et al. Adaptive image 
watermarking algorithm based on improved perceptual models [J]. AEU- 
International Journal of Electronics and Communications, 2013, 67 (8): 690- 
696. 

[12] Amin S, Hamid E. Remote sensing image matching based on adaptive 
binning SIFT descriptor [J]. IEEE Trans on Geoscience and Remote 
Sensing, 2015, 53 (10): 5283-5293. 

[13] P 'E 9], 刘 瑞 珍 . 基于 圆 形 邻 域 SIFT 特征 描述 子 的 3D 人 脸 识别 [I]. 
南京 邮电 大 学 学 报 : 自然 科学 版 , 2015, 35 (3): 11-17. (Lu Guanming, Liu 
Ruizhen. 3D face recognition based on SIFT feature descriptor in circular 

Journal 


neighborhood [J]. of Nanjing University of Posts & 


201804.02149v1 


chinaXiv 


录用 稿 


Telecommunications: NATURAL SCIENCE EDITION, 2015, 35 (3): 11- 
17.) 

[14] Fazli S, Moeini M. A robust image watermarking method based on DWT, 
DCT, and SVD using a new technique for correction of main geometric 


attacks [J]. Optik-International Journal for Light and Electron Optics, 2016, 


hinaXiv 合 作 期 和 


李 莹 莹 ， 等 : 基于 DWT 最 优 多 子 图 和 SIFT 几何 校正 的 和 鲁 梯 水 印 算 法 


127 (2): 964-972. 

[15] Li Jianzhong, Zhu Yinghui. A geometric robust image watermarking scheme 
based on DWT-SVD and Zernike moments [C]// Proc of the 3rd IEEE 
International Conference on Computer Science and Information Technology. 


2010: 367-371. 


